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ABSTRACT 

Context. The detection of small mass planets with the radial- velocity technique is now confronted with the interference of stellar noise. 
HARPS can now reach a precision below the meter-per-second, which corresponds to the amplitudes of diiferent stellar perturbations, 
such as oscillation, granulation, and activity. 

Aims. Solar spot groups induced by activity produce a radial-velocity noise of a few meter-per-second. The aim of this paper is to 
simulate this activity and calculate detection limits according to different observational strategies. 

Methods. Based on Sun observations, we reproduce the evolution of spot groups on the surface of a rotating star. We then calculate 
the radial-velocity effect induced by these spot groups as a function of time. Taking into account oscillation, granulation, activity, and 
a HARPS instrumental error of 80 cms"', we simulate the effect of different observational strategies in order to efficiently reduce all 
sources of noise. 

Results. Applying three measurements per night of 10 minutes every three days, 10 nights a month seems the best tested strategy. 
Depending on the level of activity considered, from logi^J^^ = -5 to -4.75, this strategy would allow us to find planets of 2.5 to 
3.5 Me in the habitable zone of a KIV dwarf. Using Bern's model of planetary formation, we estimate that for the same range of 
activity level, 15 to 35% of the planets between 1 and 5 Me and with a period between 100 and 200 days should be found with 
HARPS. A comparison between the performance of HARPS and ESPRESSO is also emphasized by our simulations. Using the same 
optimized strategy, ESPRESSO could find 1.3 Me planets in the habitable zone of K dwarfs. In addition, 80 % of planets with mass 
between 1 and 5 Me and with a period between 100 and 200 days could be detected. 

Key words, stars: individual: o-CenB - stars :planetary systems - stars: activity - Sun: activity - Sun: sunspots - techniques: radial- 
velocities 

1. Introduction turbation produced by three different physical phenomena: os- 

cillations, granulation, and magnetic activity. 
We have now discovered more than 400 exoplanets with the Oscillations of solar type stars, which can be seen as a di- 
radial-velocity (RV) techniquefl The majority of them ai-e very latation and contraction of external envelo pes over timescales 
massive and on very short period orbits, something that was of a f ew minutes (5 rninutes for the Sun: ISchriiver & Zwaniil 
quite unexpected from the theories of giant-planet formation. 200Q|; iBroomhall et alJ l2009l) . is caused by pressure waves 
However, since a few years, gi ant planets rn uch more similar (p-modes) that propagate at the surface of solar type stars, 
to the solar system giants (e.g. IWri^ht et all 12008) as well as The individual amplitudes of p-modes are typically from a 
super-Earth planets have been detected (mass from 2 to 10 Me) few to tens of cms"', but the interference of tens of modes 
(e.g. Mayor et al. 2009b a; Udry et al. 2007). This has become with close frequencies introduce RV vai-iations of several 
^ possible thanks to three important improvements of the RV tech- ems', de pending on the star's spe c tral type and evolution- 
■ - ' nique. First, the precision of spectrographs improved consider- ary stage (Beddins & Kieldsen"2007'; 'Bouchv & Carrier! 120031: 
ably, reaching now a 100cm s"' precision level (typically 100 Bedding & Kjeldsen 2003; Schrijver & Zwann 2000). The am- 
cms"' in 1 minute for a V = 7.5 KO dwarf, using t he HARPS plitude and period of the oscillation modes increase with mass 
spect rograph on the ESO 3.6 meter telescope, see iPe pe et al.1 along the main sequence. Theory and observations show that the 
l2005h . A second improvement has been achieved through an ap- frequencies of the p-modes rise with the square root of the mean 
propriate observational strategy, which allo ws us to average out density of the stai- and that their amplitudes are proportional to 
perturbations caused by stellar oscillations dSantos et al. 2004). the ra tio of the luminosity over the mass (Christensen-Dalsgaarl 
Finally, several years of follow-up helped us find long-period |2004tiOToole et alJI2008l) . 

planets as well as very small mass ones. At the 100 cm s"' level The convection in external layers of solar type stars drives 
of accuracy, we start to be confronted with noises caused by the different phenomena of granulation (granulation, mesogranula- 
stars themselves. Stellar noise is the result of three types of per- tion, and supergranulation), which also affect the RV measure- 
ments on time scales going from several minutes to several 

hours. In order of timescale and size of convective pattern, we 

' see The Extrasolar Planets Encyclopaedia, http://exoplanet.eu first have granulation, whose typical timescale is shorter than 25 
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minutes (iTitle et al.lll989t iDel Moroll2004 . Then comes meso- 
granulation and finally supergranulation, with timescales of up 
to 33 hours (Del Moro et al. 2004). When integrated over the 
entire stellar disk, all these convection perturbations present a 
noise level on the order of the meter-per-second. 

On longer timescales, similar to the star rotational period, 
the presence of activity related spots and plages perturbs pre- 
cise RV measurements. Spots and plages on the surface of a 
star will break the flux balance between the red-shifted and 
the blue-shifted halves of the star. As the star rotates, a spot 
group, or a plage, moves acros s the stellar disk and produces 



an apparent Doppler shift ( Saar & Donahue 
1200 it iHuelamo et all l2008t iLagrange et al. 



119971: lOueloz et al.1 



20101). This effect 



can be hard to distinguish from the signal caused by the pres- 
ence of a plan et. For the Sun at maximum activity of cycle 23, 
iMeunier et aP (12010) find a noise related to spot groups and 
plages of 42 cms"'. Because the temperature of spot groups 
and plages are different from the mean stellar surface, and be- 
cause active regions contain both, the noise induced by them 
will usually be compensated, but not entirely, because the sur- 
face ratio between s pots and plages varies (e.g. IChapman et alJ 
I2OOI.) . According to lMeunier elal] (1201 Oh . the major perturb a- 
tive effect of activity on RVs is not the one induced directly by 
spot groups and plages, but the one caused by the inhibition 
of convection in activ e regions (e.g. Dravins 1982; Livingsto^j 
ll982l:lBrandt&Solank i 1990; Gray 1992). This effect leads to a 
blueshift distortion of spectrum lines, which results in a noise 
varying from 40 cms"' at minimum activity to 140 cms"' at 
maximum. 

The Earth produces a radial-velocity perturbation of 9 cm s"' 
on the Sun, which would be completely masked by the stellar 
noise. If we manage to understand the structure of these dif- 
ferent kinds of noise, we could use appropriate observational 
strategies to reduce the stellar noise contribution as much as 
possible, and thus find very small mass planets far from their 
host star This investigation of optimized observational strate- 
gies is essential for future and more accurate instruments, such 
as ESPRESSO@VLT (http://espresso.astro.up.pt /, precision ex- 
pecte d: lOcms"') or CODEX@E-ELT (e.g. iPasquini et al] 
I2OO8I precision expected: 2 cms '), in order to detect Earth 
twins (Earth-mass planets in habitable regions). 

The two first types of noise, which are caused by oscilla- 
tions and granulation p henomena, h ave been discussed in a pre- 
vious paper (Dumusqu e et al.ll201 iL hereafter Paper I). Starting 
from HARPS asteroseismology measurements, we characterized 
these two kinds of noise using the power spectrum representa- 
tion. Then we derived an optimized observational strategy, re- 
ducing at best the stellar noise contribution while keeping a rea- 
sonable total observational time. 

In the present paper, we continue this first study and include 
the noise coming from magnetic activity spot groups. Starting 
from observations of the Sun, which is the only star where we 
can resolve spot groups, we simulate the appearance of spot 
groups on its surface and calculate the RV contribution. Adding 
the noise coming from magnetic activity to the results obtained 
in Paper I, we develop an observational strategy that simulta- 
neously reduces the three kinds of stellar noise. To finish, we 
calculate the corresponding detection limits, as well as the ex- 
pected number of planets that could be found by comparing our 
results with Bern's model of planetary formation. 



2. Simulation of different activity levels 

For decades surveys have studie d the activity level of FGK stars 
in the solar neighborhood (e.g. ' Baliunas et al.|[T995t iHall et all 
2007). One of the main results of thes e surveys is the bi-modahty 
of the activity leve l distribution (e.g. IVaughan & Prestonll 19801: 
iHenrv et"al]|1996h . About 70% of the stars, even those with a 
clear Sun-like activity cyc le, seem to have an activity level be- 
low log(/?;^^) = -4.75 (see Noyes et al. 1984, for more informa- 
tion about the chromospheric emission ratio, /?^^). The others 
have a higher activity level, and no crossing between the low- 
activity and high-activity region seems to occur. Early results of 
the Kepler mission arrive at the con clusion that th ere may be 
more active stars than expected (see Basri et al. 2010). Of more 
than 100000 stars with a surface gravity, log g, higher than 4 and 
an effective temperature between 3200 and 19000 K, 50 % seem 
to be more active than the active Sun. Nevertheless, very low- 
activity stars have been discovered in this sample. 

In the present paper, we simulate the activity of stars below 
log(R'^^) - -4.75. The goal is to reproduce as realistically as 
possible the appearance of stellar spot groups on the surface 
of these stars. Because the Sun is the only star with an activ- 
ity level between -5 and -4.75 for which we can resolve stellar 
spot groups, we will use it as a proxy for low-activity stars. Stars 
with a logiR'fjf,) higher than -4.75 will not be simulated. First 
of all because of a lack of information about the activity phe- 
nomenon of these stars, and secondly, because the activity noise 
will probably be too important to find very small mass planets in 
habitable regions with the RV technique. 



2.1. Activity of the Sun 

The Sun has a 22-year magnetic cycle as well as a 11 -year 
sunspot cycle. During these 11 years, the activity level varies 
between log(/?^^) - -5 at minimum and log(/?^^) - -4.75 at 
maximum. At minimum there is no spot group on the surface of 
the Sun. When the activity level rises, several spot groups start to 
appear at a latitude of about +30°. Before reaching the maximum 
activity level, the number of spot groups will increase progres- 
sively, whereas the latitude of spot groups will decrease (+15° 
at maximum activity). After this maximum phase, the migration 
of spot groups toward the equator continues and the number of 
groups decreases progressively until they all disappear At the 
end of a cycle, the activity level returns to log{R'^j^) = -5, and 
another sunspot cycle of 1 1 -year begins again. 

Because the latitude of the spot groups changes during the 
sunspot cycle, it is important to recall that the Sun has a differ- 
ential rotation in latitude. Therefore, spot groups rotate faster at 
the equator than at the poles, according to the equation 



cj - A -i- B sin^ 6, 



(1) 



where co is the angular speed in degree /day, 6 the latitude, and 
A and B a re two constan ts equal to 14.476 + 0.006 and -2.875 ± 
0.058 (seeE owar dl[T996l) . 

During the 11 -year sunspot cycle or an even longer time, 
the spot groups do not appear randomly in longitude, but 
l^deg around preferred activ e longitudes (e.g. Ilvanovl 120071: 
i Berdvugina & UsoskinI 120031) . These preferred longitudes turn 
with the Sun's differential rotation and are observed for large 
spot groupfl Observation of the Sun over several sunspot cycle 
has shown that 



^ Large spot groups represent 11% of the total number of spot groups 
observed, corresponding to 67% of the total spot group area. 
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- there are normally two active longitudes in the northern 
hemisphere and two in the southern one, 

- the two active longitudes are shifted by 180''*^* in each hemi- 
sphere, 

- the active longitudes in the south are shifted by 90''*"* in re- 
gard to those in the north, 

- sometimes, each hemisphere can have a third active longi- 
tude, which is shifted by 90'''°* with regard to the two normal 
ones. 

When sunspots appear, their diameter is very small, approxi- 
mately 0.003 Rq. The majority of them will disappear in a couple 
of hours, but a few will evolve and become bigger spots, which 
will gather together as spot groups. These ca n last for severa l 
days to a few months. According to Howard & MurdinI (l2000h . 
the lifetime distribution of spot groups is 

- 50% of the spot groups have a Ufetime between 1 hour and 2 
days, 

- 40% of the spot groups have a lifetime between 2 days and 
1 1 days, 

- 10% of the spot groups have a lifetime between 1 1 days and 
2 months. 

- For the Sun at maximum activity (log(/?^^) - -4.75), there 
is one group spot per sunspot cycle that can last from 2 to 6 
months. 

During its lifetime, a spot group will rapidly increase in area un- 
til the maximum size is rea ched, after which the decr ease will 
follow in a slower way (e.g lHoward & MurdinI l20()0h . The in- 
crease and decrease times are not precisely defined. For our sim- 
ulation, we will set these times to 1/3 and 2/3 of the total lifetime 
of a spot group. 

The maximum spot group filling factofl fs.max, is linked 
to the lifetime of the spot group, T, by the equation 
(iHoward & Murdinll2000l) 



fs,n 



(2) 



where fs,max is expressed in hemisphere and T in days. 

Once the lifetime and the maximal size of spot groups are 
defined, we need a probabilistic process to mimic the observed 
appearance of spot groups on the solar surface. According to 
[Hoyt & Schatten ( 1998), the appearance of spot groups can be 
described by a Poisson law. We can therefore define the prob- 
ability, P, to have a certain number of spot groups, A^, by the 
formula 



-Xt 



pmt + T)-N{t))^k] = 



A: = 0, 1,- 



(3) 



where f is the time, t the time step and A the average appearance 
of spot groups per unit of time. 

The upper panel in Fig. [1] represents the continuous 
daily observation of sunspots (Marshall Space Flight Centre, 
http://solarscience.msfc.nasa.gov/SunspotCycle.shtml). Clearly, 
the number of spots varies with the 11 -year sunspot cycle of 
the Sun. Comparing the two graphs in Fig. [1] we notice that 
the number of spots is correlated with the S -index and thus 
wit h the \o^( R'ify) (because S and log(/?^^) are proportional, 
see iNoves et al. (1984)). In conclusion, no spot corresponds to 
\og(R'fjj^) - -5, and the maximum number of spots corresponds 
to \og(R'jjj^) - -4.75. According to this argument, we decide to 
simulate three diff'erent levels of activity: 
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^ The filling factor is just the ratio of the area of the spot group over 
the area of the visible hemisphere. 



Fig. 1. Upper Panel: Monthly averages of the sunspot 
numbers. We clearly see the 11 -year sunspot cycle (from 
http://solarscience.msfc.nasa.gov/SunspotCycle.shtml). Lower 
Panel : Variation of the S-ind ex as a function of time for the Sun 
(from lLockwood et al.ll2007h . 



- log(/?^^) - -5: the Sun is at minimum activity and there is 
no spot group on the solar surface. Therefore, no noise from 
the spot groups perturbs the RVs (case of Paper I). 

- \ci%(R'ff^ - -4.75: the Sun is at maximum activity. The 
number of spots in the visible hemisphere is approximately 
equal to 150 and they a re located at +15° of latitude (e.g 
IHoward & Murdinll2000l) . 

- log(7J^^) = -4.9: the Sun is at medium activity, which cor- 
responds to a total number of 75 spots at +22.5° of latitude 
in the visible hemisphere. 

The activity level is given by a defined number of spots, but 
the properties defined above, such as the lifetime distribution, 
the maximum filling factor (Eq. |2]l, and the appearance statis- 
tical law (Eq. |3]l, are given for spot groups. In order to use all 
these properties, we need to know the number of spot groups 
of a defined number of spots. This can be done using the work 
of Hathawav & Choudharv (2008), which gives the number of 
sunspots per spot group as a function of the spot group maxi- 
mum filling factor (see Fig.|2]i. 



2.2. Activity simulation 

In order to check our ability to detect Earth-like planets through 
activity effects, our simulations will try to mimic the Sun activity 
phenomena as much as possible. 

For the simulation, we start with the following initial condi- 
tions: 

- We set the number of spots, A^,, for a given activity level. 

- Then, according to the lifetime distribution described in Sect. 
12.1! we create a first spot group and calculate its maximum 
filling factor using its lifetime (Eq.|2]i. 
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Fig. 2. Number of sunspots per spot group (right scale) as a func- 
tion of t he spot group maximal area in millionth of hemisphere 
(source: iHathawav & Choudharvl2008l) . 





time [d] 



Fig. 3. RVs produced by a spot group at 15"'^^ of latitude, with a 
filling factor of 0.001 . The simulation is made using the program 
SOAP. 



- This maximum filling factor gives us, using Fig. |2l the 
number of spots present in this spot group, Ns, group- If 
Ns, group < Ns, we create a second spot group and so on, until 

^YjNs, group ~ Ns- 

- Each spot group is put on the latitude corresponding to the 
activity level and on the different active longitudes, with no 
possible overlap. The longitudes of the spot groups are se- 
lected using a Gaussian distributiorQ centered on active lon- 
gitudes and with a dispersion of 15'''°*. 

Once the initial conditions are created, we start the time evo- 
lution of the system, using a rotational period of the Sun at the 
equator of 26 days. The time step of the simulation is set to two 
hours. All spot groups begin with a null area and increase their 
size linearly until they reach the maximum after 1/3 of their life- 
time (see Eq.|2]). The decreasing phase is linearly during the re- 
maining 2/3 of their lifetime. The number of new spot groups at 
each time step is given by the appearance statistical law (see Eq. 

At the beginning of the simulation, the value for the aver- 
age appearance of spot groups. A, is unknown. To find the right 
value, we begin the simulation with an arbitrary number and cal- 
culate the total filling factor for each time step. Using Fig.|2] this 
filling factor give us the total number of spots. If the mean total 
number of spots in not equal to A^^ in several years, we restart the 
simulation with another value for A. Finally, after several tests, A 
is set to 2.4 and 4.8 for an activity level of log(/?^^) = -4.9 and 
\og{R'fjj^) - -4.75, respectively. 

For each level of activity, we generate 100 simulations with 
an equal probability of having two or three active longitudes per 
hemisphere. 



2.3. Radial-velocities induced by activity-related spot groups 

When spot groups are present on a rotating star, they break the 
flux balance between the red-shifted and the blue-shifted halves 
of the star As the star rotates, a spot group moves across the 
stellar disk and produces apparent Doppler shifts. This modula- 
tion is often hard (sometimes impossible) to distinguish from the 
Doppler modulation caused by the gravitational pull of a planet. 




Because 68 % of the values generated by a Gaussian distribution 
will fall in the range yu ± cr, where fi is the mean and cr the dispersion, 
this technique allows us to put 68% of the total spot groups area at 15° 
around the active longitudes (67 % observed, see Sect. 12. It . 



0.0000 0.0002 0.0004 0.0006 0.0008 0.0010 0.0012 0.0014 0.0016 

Fig. 4. Semi-amplitude of the RV signal, Krv induced by a spot 
group as a function of the spot group filling factor These results 
have been obtained with the program SOAR We can see in black 
the result for spot groups at a latitude of 15"'^*, and in gray (red), 
of22.5''''«. 



To model the radial-velocity variations imposed by stellar spot 
groups, we used the program SOAP (Bonfils & Santos, in prep.). 
SOAP computes the rotational broadening by sampling the stel- 
lar disk on a grid. Each grid cell is assigned a Gaussian func- 
tion that represents the typical emerging spectral line (or, equiv- 
alently, the spectrum's Cross Correlated Function (CCF)). All 
cells are Doppler-shifted according to their projected velocities 
toward the observer's line-of-sight, and averaged with a weight 
following a linear limb-darkening law (a - 0.6). When spot 
groups are added to the stellar surface, SOAP computes which 
grid cells they occupy and change the weight for those given 
cells. The weight can be set to zero for a dark spot group, to a 
fraction of the stellar brightness if we assume a given brilliance 
for the spot group, or to more than the stellar brightness to simu- 
late plages. Finally, for a given spot group configuration, SOAP 
delivers the modified spectral line, the flux, the radial-velocity, 
and the bisector span as a function of the star's rotating phase. 
Note that SOAP merely calculates the flux effect induced by spot 
groups. It does not model any type of granulation, therefore, the 
inhibition of convection (granulation) in active regions s uch as 



spot groups is not took into account (iMeunier et 



regions su( 
al] l20I0h . 
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Fig. 5. Simulations for an equatorial rotational period of 26 days and an activity level log(/?^^) - -4.9 (left panel) and for 
\og{R'^j^) - -4.75 (right panel). The horizontal line on the top and bottom graphs are the RV mean and the spot number mean 
over the 4 years of the simulation, respectively. Just the 400 first days are shown here for clearness. 



Normally a spot group is constituted of several small spots. 
For the simulation we considered each spot group as a unique 
big spot to facilitates the operation. We used the program SOAP 
to calculate the RV effect induced by a spot group presenting 
a fining factor, /,i of 0.001 and a latitude of 15''''« (Fig. [3]) or 
22.5''"^. The first latitude corresponds to a \og{R'jjf,) of -4.75 and 
the second to a log(/?^^) of -4.9. To get the RV effect for any spot 
group filling factor, fs2, we first calculate the RV semi-amplitude 
with SOAP, Krv, which is induced by different spot group sizes. 
Then we fitted a power law model to these data (Fig. HI and 
found that 

- for a spot group at a latitude of 15''''*, Krv °^ fs '^^, 

- for a spot group at a latitude of 22.5''''*, Krv °^ fs '^'^- 

These results ar e con sistent with the relation presented in 
ISaar & Donahud (Il997i) . where the authors have found that 
Krv y?'^"- We will use this relation in our simulation for 
all latitudes. Supposing that besides the amplitude, the RV sig- 
nal has the same shape for any spot group size, we calculate 
the proportionality coefficient between two signals, which is 
Krvi/Krvi - ifsilfsif'^- Multiplying the RV effect induced by 
a spot group with - 0.001 by this coefficient yields the RV 
effect induced by any spot group size. 

The simulation does not consider spot groups with a lifetime 
shorter than one day, because the maximum filling factors will 
be very small an therefore the RV contribution negligible. 

We present here a simulation for the Sun, which is seen equa- 
tor on, therefore the inclination angle between the line of sight 
and the rotation axis of the star, /, equals 90°. For other stars, this 
assumption cannot be made because for the majority of them we 
do not have access to the / angle. Moreover, when / decreases, 
the RV which is equal to v sin / decreases as well . Thus the RV 
effect of spot groups crossing the visible stellar disc, which is 
proportional to v sin /, will be lower. In conclusion, considering 
i = 90° in the simulation gives us the maximum RV effect that 
spot groups can induce. 

2.4. Comparison with other works 

An example of simulations for log(7?^^)=-4.9 and -4.75 can be 
seen in Fig. |5] The amplitudes found can be compared to a re- 
cent result from Meunieretal. (2010), who calculated the RV 
induced by the real distribution of spot groups during the solar 
sunspot cycle 23. For maximum activity, they find a RV variation 



of 48 cm s"' . Note that this variation is valid only for spot groups 
and was derived without taking into account plages or inhibition 
of the convection in active regions. In our case, after 100 simu- 
lations of a high-activity level, we arrive at a valu e of 63 cm s~'. 
Our estimate is 30 % higher that the one found by'Meunier et ajj 
(2010), which is based on real position and size of spot groups. 
This can be explained because iMeunier et al.l (12010 ) used ob- 
servations of the solar cycle 23. This was a fairly quiet cycle, 
presenting at maximum no more than 100 sunspots in the visi- 
ble hemisphere (see Fig. [T]). In our case, the maximum number 
of sunspots was based on cycles 21 and 22, which were more in- 
tense (150 sunspots in the visible hemisphere at maximum activ- 
ity level). In Fig.|7]we show the calculated RV variation obtained 
for different values of the total spot number A quadratic fit to the 
data points indicates a RV variation of 5 1 cm s ~' for 100 spots , 
which is very close to the value cited in Meunier et al.l (120101). 
Thus, our simulation, starting from very simple characteristics 
of spot groups such as the lifetime distribution, the spot number 
according to the activity level, and the Poisson process of spot 
group appearance, manages very well to represent the RV effect 
induced by spot groups. Note that the inhibition of convection 
in active r egions is not yet incl uded in our simulation, and ac- 
cording to iMeunier et al.l (|2010'), this effect could be dominant. 
Therefore, the RV effect of the activity is not yet fully simulated 
and could be underestimated. The very good agreement compar- 
ing the spot flux effect between real Sun observations and the 
simulation indicates that our simple model of spot groups is re- 
alistic. Including the inhibition of convection in the simulation is 
in progress and will be the topic of a forthcoming paper. Our ac- 
tivity si mulation presents an a dvantage compared with the work 
done by iMeunier et al.l (1201 Oh in the sense that we can adapt it 
to other stars than the Sun, if basic characteristics of the spot 
groups are known. 

According to lOueloz et af] (|2009|) and iBoisse et al] (1201 Ol) . 
the RV signal caused by stellar spots is found at the rotational 
period of the star, Proi, and its two first harmonics. When con- 
sidering few big spots at the surface of the star that can change in 
size but not disappear rapidly, this assumption is valid because 
the activity signal stays in phase with rotation. In the two pa- 
pers cited above, the authors consider active stars that are known 
to be dominated by few big spots that can live for several ro- 
tational periods. In our simulation, were we model non-active 
stars, we consider a lot of spots that are allowed to disappear 
rapidly and reappear at other places any time. This random be- 
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spot number 

Fig. 7. Evolution of the RV variation as a function of the total 
number of spots. Each point corresponds to 100 simulations. 



havior, ruled by a Poisson law, kills the activity signal phase even 
on timescales similar to the rotational period of the star. The RV 
signal however keeps a kind of periodicity on a rotational pe- 
riod timescale, but the period can be different form Proi and its 
two first harmonics. In Fig. |6]we plot the periodograms and the 
RVs for 31 day slices (1.2 x P^,) of the simulation showed in 
Fig. |5] with an activity level of log{R'^^) = -4.75. Looking at 
the periodograms, we see that sometimes important peaks are 
pr esent at Pmt, Prnt/2 a nd P ,„t/3 (plot on the rig ht), as suspected 
bv lOueloz et akl (120091) and lBoisse et al. At other times, 

the peak at Proi is clear but nothing shows up at P,o/ /2 and Proil'i 
(plot on the left) or at the rotational period and the two first har- 
monics (plot on the center). Looking now at the RV plots, we see 
that the model of three sine waves with period P/„, Pfii/2 and 
Pfitl^ can fit the simulated RVs nicely. The comparison between 
the rms of the simulated RVs and the corrected one (simulated 
RVs minus fitted sine waves) shows that the three sine waves 
model reduces the activity noise by 70 %. Therefore this model 
can be used to correct short-term activity noise. However, if a 
small mass planet orbits the stars with a period close to 
the three sine waves model will kill this signal too. We have 
therefore to be very careful when using these models, and only 
the study of other CCF parameters such as the full width at half 
maximum (FWHM) or the bissector span (BIS) will be able to 
give us clues on the nature of the signal: short-term activity or 
planet. This work will be done once the inhibition of the con- 
vective blueshift in active regions will be implemented into our 
model. In addition, the estimate of th e rotational period us ing the 
three sine waves model, as shown in lOueIoz et al.l (l2009l) . is no 
longer valid. Indeed, we see that the fitted period, f /,-,, can vary 
between 22 to 34 days, whereas the true rotational period of the 
star is fixed to 26.4 days in the simulation. 

3. Efficient observational strategies 

The RV signal induced by a spot group can be described as the 
succession of null contribution when the spot group is in the hid- 
den hemisphere, and as sinusoidal shapes when the spot group 
moves across the visible hemisphere (Fig. 13). This RV signal is 
regular in time and is composed of sinusoidal signals with peri- 
ods egualtoflierotational period of the star and its first harmon- 
ics (IBoisse et al I2OIOI) . If there are several spot groups, these RV 
signals will add or cancel themselves, leading to a quasi-periodic 
signal with unpredictable but limited amplitude. Because the pe- 
riod of the RV signal varies between approximately 10 to 30 



days, the idea is to average out these efi'ects by an adequate sam- 
pling of the variation time scale. 

In Paper I, a good observational strategy to average out os- 
cillation and granulation noise was to take three 10-minute mea- 
surements per night with two hours of spacing. This strategy was 
applied to the actual calendar of HD69830, which is one of the 
most followed stars by HARPS. This calendar, besides obser- 
vation gaps owing to bad weather, is very similar to 10 consecu- 
tive nights of observation every month. This idealistic strategy of 
three measurements per night of 10 minutes each, on 10 consec- 
utive nights every month, will be referred to below as the 3N1 
strategy. Because the typical time scales of activity noises are 
longer than 10 days (see Fig. |5]l, the 3N1 strategy is not opti- 
mized to average out this kind of perturbation. In order to aver- 
age out all three kinds of noise at the same time, we have to cover 
a period of a month with measurements, while keeping the ben- 
efit of the higher frequency sampling. Taking into account this 
last argument, we simulate the effect of two new observational 
strategies. Keeping the three measurements per night of 10 min- 
utes and the 10 nights of observation per month, we change the 
separation between each observational night. The new consid- 
ered observational strategy will measure the star every second 
night (hereafter 3N2 strategy) and every third night (hereafter 
3N3 strategy). For the 3N1, 3N2, and 3N3 strategies the calen- 
dar of HD69830 was not suitable any more and we had to create 
new calendars, asking for the following criteria: 

- The total observational time span is four years. 

- The star disappears from the sky during four months each 
year 

- Each month of observation, the beginning of the measure- 
ments is randomly set. 

In addition, we randomly suppress 20 % of nights to take into 
account bad weather or technical problems. These calendars rep- 
resent a total of 256 nights of measurements over a time span of 
four years. In Fig. |8] (left panel) we can compare these diff'er- 
ent observational strategies. The RV variation, induced by spot 
groups for a maximum level of activity (log(7?]y^) = -4.75), 
reaches a level of approximately 60cms"^ By using the 3N3 
strategy, we manage to reduce this level of RV variation to 
10 cm s ' when averaging the effect in 10 day bins. For the same 
binning, the 3N1 strategy reduces this level to 30 cms"'. 

We can now calculate nnsnvb as a function of the binning in 
days, including all types of noise: oscillation, granulation, activ- 
ity, and instrumental (80cms"'). The result for aCen B (KIV) 
can be seen in Fig.[8](right panel). 

We notice that the 3N3 strategy averages out much more ac- 
tivity noise after binning than the 3N1 strategy. Without activity 
(logiR'^f,) = -5) the two strategies are identical. The RV varia- 
tions reachable using 10 days bins (suitable for long period plan- 
ets) with the best considered strategy are 22 cms"', 25 cms"', 
and 28 cm s"' for a log(/?^^) equal to -5, -4.9 and -4.75, respec- 
tively. If there is high activity, the improvement brought by this 
optimal strategy reaches 30 % compared to the 3N1 strategy. 

In conclusion, sampling the complete rotational period of the 
star with the same number of measurements, but with more sep- 
aration between them, averages out the activity noise to a con- 
siderable extent. Thus, the strategy needs to be adjusted to the 
stellar rotational period. Here for the Sun, which has a 26 day 
rotational period, the 3N3 strategy is the optimal one with 10 
measurements per month. For a star with a 20 day rotational pe- 
riod, the 3N2 strategy would average out the stellar noise better 
In addition, increasing the gap between observational nights by 
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Fig. 6. Periodograms (top) and corresponding RVs (bottom) for three different slices of 3 1 days for the activity simulation shown in 
Fig. |5](log(/?^j^) = -4.75). In the periodograms the horizontal line corresponds to the 1 % false-alarm probability and the vertical 
dashed lines to the rotational period of the star, P,o,, as well as P,of/2 and P,of/3. On the RV plots, the dots correspond to the 
simulated RVs and the triangles to a fitted function with three sin waves of period P/,-,, Pfii/2 and P fit/3 (Pfu is a free parameter). 
Prot as well as the rms of the simulated RVs can be found in the upper left part of the plot. Py„ as well as the rms of the corrected 
RVs (simulated minus fitted sine waves) can be found in the upper right part of the plot. 




Fig. 8. Evolution of the binned RV variation, rmsRvb, owing to activity-related group spots (left panel) and to all kind of stellar 
noises (right panel) as a function of the binning in days. The three strategies shown for different log(R'^^) (given for each lines) use 
three measurements per night of 10 minutes and 10 nights of observation per month. The only difference lies in the spacing between 
each night of observation. Thicks lines correspond to the 3N1 strategy, thin lines to the 3N2 strategy, and dashed lines to the 3N3 
strategy. 



more than three days for slow rotators is not recommended, be- 
cause it will induce a poor sampling of short period planets. 

4. Detection limits in the mass-period diagram 

In order to derive detection limits in terms of planet mass and pe- 
riod accessible with the different studied strategies, we calculate 
the false-alarm probability (FA P) of detection usi ng bootstrap 
randomization (Endl et al. 20"0H lEfron & Tibshiran i' 1998). We 
first simulate a synthetic RV set containing stellar, instrumen- 



tal, and photon noises as we did in the previous sections. The 
calendars used depend on the selected strategy, but all have a 
total of 256 observational nights over a time span of four years 
(see Sect.O. Then we carry out 1000 bootstrap randomizations 
of these RVs and calculate the corresponding periodograms. For 
each periodogram, we select the highest peak and construct a 
distribution of these 1000 highest peaks. The 1 % FAP corre- 
sponds to the power, which is only reached 1 % of the time. 
The second step consists in adding a sinusoidal signal with a 
given period to the synthetic RVs. We calculate the periodogram 
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Fig. 9. Detection limits for aCen B(KIV) for different activity levels and strategies. The thick line at the top corresponds to the 
present HARPS observational strategy (one measurement per night of 15 minutes, 10 consecutive days each month) for an activity 
level log(/?^^) - -4.75. The second thick line at the very bottom corresponds to the detection limit, using ESPRESSO, for the 3N3 
strategy and an activity level log(/?^^) = -5. Besides these two thick lines, from top to bottom, the dashed thin lines correspond to 
the 3N1 strategy using HARPS, for an activity level log(7?^^) - -4.75, -4.9 and -5, respectively. The 3N3 strategy using HARPS 
corresponds to the continuous thin lines on the same order of activity level. Finally, the small dots represent the planets expected 
from Bern's model with a random sin / for each body. We see that the cut-off in period around 12 days of the model introduces an 
anomaly in the planet distribution. The habitable zone is derived using the model of [S elsis et al. (2007) applied to the temperature 
and the luminosity of a Kl V dwarf, calculated using the results of lde Jager & Nieuwenhuiizeni(il987i) . 



of these new RVs and compare the height of the observed peak 
with the 1 % FAP We then adjust the semi-amplitude of the sig- 
nal until the power of the peak is equal to the 1 % FAP. The 
obtained semi-amplitude corresponds to the detection limit of a 
null eccentricity planet with a confidence level of 99 %. To be 
conservative, we test 10 different phases and select the highest 
semi-amplitude value. 

In practice, we simulate 100 RV sets for each activity level 
(log(/?^j^)=-5, -4.9 and -4.75) and strategy, and to be conser- 
vative, only the 10 "worst" cases per strategy were considered. 
It is important to note that for each strategy we simulate the 
RVs according to the respective calendar (10 measurements per 
months on consecutive nights or every second or third night, 
eight months a year over four years) and we removed randomly 
20 % of the nights to simulate bad weather or technical prob- 
lems. The observational calendars are therefore realistic. Figure 
|9] presents for different levels of activity the 99 % confidence 



level detection limits expected for the 3N1 and the 3N3 strate- 
gies using HARPS. 

The thick line in top part of Fig. |9] corresponds to the de- 
tection limit for a maximum-activity level (log(/?^^) = -4.75) 
using the present HARPS observational strategy (INl strategy 
= one measurement of 15 minutes on 10 consecutive days per 
month, see Paper I for details). Obviously, there is a great im- 
provement brought by the 3N1 and 3N3 strategies compared to 
what is done at present on HARPS. 

We can also compare the level of the detection limits ob- 
tained for the 3N1 and 3N3 strategies. For example, at a 100 
day period, the improvement brought by the 3N3 strategy is 
4%, 10% and 17%, for an activity level of log(/;^^) = -5, 
-4.9 and -4.75, respectively. As expected, the improvement in- 
creases with a rising activity level because the 3N3 strategy is 
optimized to reduce activity noise. For minimum activity, which 
corresponds to no spot groups and consequently no noise from 
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activity in our simulation, the detection limits obtained with the 
two strategies should be identical. The 4 % difference can be ex- 
plained by a more regular sampling (over all months with the 
3N3 strategy and only 10 consecutive days per month for the 
3N1). In conclusion, the 3N3 strategy seems to be an efficient 
strategy to average out all kinds of noise. For example, for a 
period of 200 days, which corresponds to the habitable region 
around a Kl dwarf such as a Cen B , the detection limits for the 
3N3 strategy ranges between 2.5 and 3.5 depending on the 
activity level. Following these results, it seems clear that even in- 
cluding the activity noise related to spot groups, planets smaller 
than 5 in habitable regions could be detected with HARPS 
using an appropriate observational strategy. 

Future instruments, like ESPRESSO at the VLT, will reach 
a precision level of 10 cms"'. The thick line at the very bottom 
of Fig. |9] shows the improvement when we change the instru- 
mental noise from 80 to 10 cms"'. Evidently the improvement 
is significant, reducing approximately the mass detection limit 
by 1 M®. We note, however, that because we used HARPS data 
as a basis for our simulations, it is possible that we overestimate 
the noise of the ESPRESSO-like radial- velocities. These latter 
values must then be seen as upper limits. 

5. Planet detection according to Bern's model 

The detection limits give us an idea of the smallest planet that 
could be found with an appropriate strategy. Nevertheless, this 
does not inform us on the proportion of small mass planets that 
could be detected. This requires a model of planetary formation 
that can predict the "true" population of sm all planets. We used 
Bern's model (e.g. iMordasini et al.ll2009allbl) . which is a good 
proxy of planet population. In their simulation they solve as in 
classical core accretion models the internal structure equation 
for the forming giant planet, but include at the same time disk 
evolution and type I and II planetary migration. 

Using the detection limits derived in Sect. H) we calculate 
the number of planets predicted by Bern's model that are above 
a given detectable level. The proportion of expected planets be- 
tween 1 and 5 Me that could be found with the 3N3 strategy on 
HARPS is represented in the left panel of Fig. [10] In this panel, 
we clearly see that the proportion of detection decreases when 
the activity level increases. For the period range between 100 and 
200 days, the proportion of planets with IMg < Me sin / < 5Me 
that would be found with HARPS changes from 35 % to 15 % 
when the log(/?^^) varies from -5 to -4.75. Although the ac- 
tivity greatly reduces the number of detectable planets, it is still 
possible to find some in the habitable region of K dwarfs (200 
days) for the highest case of activity (log(7?^^) - -4.75 for low- 
activity stars, see Sect.|2]l. 

We note for a similar activity level that the value of the bin 
for 300 to 400 days is lower than the one for 400 to 500 days. 
Because the mass detection limits increase with period, we ex- 
pect the opposite. This effect is owing to the observational cal- 
endar used. Indeed, because stars in the sky can generally be 
followed during eight months per year before they disappear (an 
effect taken into account in our observational calendar), the one 
year period is not well sampled, which complicates the detec- 
tion. This is well illustrated in Fig. |9] were a peak near one year 
appears for each detection limit. 

A comparison with what could be found using ESPRESSO 
is shown in the right panel of Fig. [10] If the activity level is set 
to log(/?^^) = -5, ESPRESSO could find 80 % of the planets in 
a mass range from 1 to 5 Me and a period between 100 and 200 
days, whereas HARPS could only find 35 % of them. 



6. Detection of a simulated planet 

To show that the calculated detection Umits are realistic, we 
present in this section an example of planet detection. For 
a Cen B(Kl V), we first create some synthetic RVs including all 
types of noise, oscillations, granulation phenomena, and activ- 
ity. The strategy used is the 3N3 and the activity level is set to 
log(7?^^) - -4.9. Secondly, we add the RV signal induced by a 
small planet in the habitable zone (M sin / = 2.5 Me, a 208 day 
period). To check whether the planet is detected or not, we calcu- 
late the periodogram and the 1 % and 0.1% FAR using bootstrap 
randomization (see Sect.|4]l. 

In Fig. [TT] top panels, we show the raw RVs, stellar noise 
plus planet, for the 3N3 strategy as well as the added planet. 
The periodogram with the 1 % FAP for each period is also repre- 
sented. The peak at 200 days is far above the 0.1% FAP, which 
proves that the planet is easily detected. To compare this with 
the actual measurements made on HARPS for the high-precision 
program, we apply the same process for the INl strategy (IN 
strategy in Paper I). This strategy consists of measuring the star 
15 minutes per night on 10 consecutive days per month. The bot- 
tom panels of Fig. [TT] show the results. Evidently here the planet 
is far from being detectable. 

Looking at Fig. [9] we see that the mass detection limit for 
208 days of period and for an activity level log(/?^^) = -4.9 is 
more than 3 M® for the 3N3 strategy. In this example, we see 
that using the same strategy and configuration, a 2.5 earth mass 
planet could easily be detected. Thus the detection limits derived 
in Sec.]4]are very conservative. This is because we consider only 
the 10 "worst" cases when deriving the detection limits. 

7. Concluding remarks 

In Paper I we had proposed an efficient observational strategy to 
reduce the stellar noise generated by oscillation and granulation 
as much as possible. This strategy, requiring three measurements 
per night of 10 minutes over 10 consecutive days each month, 
improves the averaging out of stellar noise coming from oscil- 
lation and granulation by 30 %, with an observational cost only 
multiplied by a factor 2. In the present paper, we add the noise 
induced by activity-related stellar spot groups. 

To study the radial-velocity (RV) effect caused by stellar 
spot groups, we consider three different activity levels. Based 
on Sun observations, we simulate the minimum solar activity 
level (logiR'fjf,) = -5), the maximum one (log(/?^^) = -4.75), 
and an intermediate one (log(/?^^) - -4.9). Comparing the RV 
variation at maximum activity given by our simulation with the 
one calculated by real position and size of spot groups from cy- 
cle 23 (Me unier et al.. ,2010). we find a very good agreement, 
51 cms"' and 48 cms"', respectively. We note that the RV ef- 
fect of activity is not fully simulated because the inhibition of 
con vection in active reg ions is not yet implemented. According 
to Meunier et al.l (1201 Oh . this effect could be important and con- 
sequently, the detection limits calcul ated here could be u nder- 
estimated. Compared to the results of iMeunier et al ] (120 1 Oh . our 
simulation is more general in the sense that we can predict the 
RV effect of activity related spots groups for other stars, know- 
ing a few characteristics (mean spot number, distribution of spot 
groups lifetime, presence of active longitudes, differential rota- 
tion). This will be very important when a better knowledge of 
activity phenomena of other stars than the Sun will be known. 
The Kepler mission should give us some clues in the near future. 

Modeling the short-term activity with a three sine waves 
function with period P/„, P/h/2 and Pfii/3 can greatly reduce 
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Fig. 10. Left panel: Proportion of expected planets with 1 < M sin / < 5 that could be found with HARPS around early-K 
dwarfs, using the 3N3 strategy. For each bin from top to bottom we have the detection proportion for log(/?^^) equal to -5, -4.9 
and -4.75. Right panel: Proportion of expected planets with 1 < Msin / < 5 that could be found with HARPS (dark or blue) 
and ESPRESSO (gray or cyan) around early-K dwarfs, using the 3N3 strategy. The activity level is set to a log(7?^^) equal to -5. 



0.12 - 




Fig. 11. Detection simulation for a 2.5 M© planet in the habitable zone {P - 208 days) of a Kl dwarf presenting an activity level 
\og{R'fjj.) - -4.9. Top panels: Detection simulation for the 3N3 strategy. On the left, we show the raw RVs including the stellar 
noise (oscillation, granulation phenomena, and activity) plus the injected planet (red curve). In the middle, the periodogram of the 
stellar noise only (red filled curve) and the stellar noise plus planet (black line) is represented. The horizontal lines correspond from 
top to bottom to the 0, 1 % and the 1 % FAP Finally we show on the right the raw RVs binned on 1 month and folded in phase with 
the period of the planet. Bottom panels: The same but for the INl strategy (strategy presently used for the high-precision HARPS 
program, one measurement per night of 15 minutes on 10 consecutive days per month). 



the spot-induced noise by approximately 70%. However, we 
have to be very careful with this method because P/„ can vary 
in a non negligible way from the rotational period of the star. 
Thus, signal of small mass planets with period similar to the ro- 
tational period of the star will be killed by this type of model. 
Only a study of other CCF parameters such as the FWHM or the 
BIS can give us clues on the true nature of the signal: short-term 
activity or planet. 



After simulating the effect of several observational strate- 
gies, the most efficient one to average out all kinds of noise is 
the 3N3. This strategy consists in measuring the star three times 
10 minutes per night, with a spacing of two hours. Then the mea- 
surements are taken every third night, 10 days every month. With 
this strategy it would be possible to find planets of 2.5 to 3.5 
Earth mass with HARPS in the habitable region of K dwarfs 
(200 days). The first mass value corresponds to a case without 



Dumusque et al.: Activity noise and planetary detection 



11 



activity (log(7?^^) - -5) and the second one to a maximum ac- 
tivity level (log(R'^fr) = -4.75). Even if the activity caused by 
spot groups introduces a non negligible noise, small mass plan- 
ets in habitable regions could be detected with HARPS with an 
appropriate observational strategy. 

Using the population of low mass planets predicted by Bern's 
model, we calculate the proportion of planets that could be found 
with the 3N3 strategy. Using HARPS, it would be possible to 
find 35 % of the planets below 5 with a period between 100 
and 200 days and a low-activity level (log(7?^^) = -5). This 
value decreases to 15 % for a log(7?^^) equal to -4.75. If we 
trust this model of formation, which is the most compatible with 
observations, HARPS could discover several planets below 5 
in the habitable region of early-K dwarfs. 

In our simulation, the 3N3 strategy appears to be the best one 
to average out activity noise. This is because of the rotational 
period of the Sun, which is fixed in our simulation at 26 days. 
Indeed, when choosing 10 observational nights a month, the best 
way to sample the entire rotational period is to observe the star 
every third night. For shorter rotational periods, increasing the 
measurement frequency would lead to a better averaging. For 
longer rotational periods, reducing the measurement frequency 
is not recommended because short period planets will be poorly 
sampled. 

Taking three measurements of 10 minutes per night gives 
us a very low photon noise when binning the data over the 
night. Thus, the photon-noise is not a limitation any longer, 
and 2 to 4-meters class telescopes can still be used to go down 
to the 10 cms ' level on bright stars. However, spectrographs 
must also intrinsically reach this level of precision and for the 
moment, only ESPRESSO© VLT (http://espresso.astro.up.pt/) is 
designed to reach this goal. With such a precision, improve- 
ments would be remarkable. The detection limits would go down 
in mass by a level of 1 M®, reaching 1.3 in the habitable 
region of early-K dwarfs. In addition, ESPRESSO would find 
80 % of the planets between 1 and 5 M@ and with a period be- 
tween 100 and 200 days, which is more than twice as much 
as w hat HARPS could find. Moreover, CODEX@E-ELT (e.g 
iPasquini et al . 2008), would improve the long term instrumen- 
tal precision down to a few cms ' level and would detect 1 
planets in habitable regions. 
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